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In addition to its global South-North anisotropy ? , there are two other major seismological 
observations relative to the Earth's inner core: asymmetry between the eastern and west- 
ern hemispheres ? ' ? ? ' ? ? and a layer of reduced seismic velocity at the base of the outer core 
?,?,?,?,?,?,?_ jjjj s 250 km thick layer has been interpreted as a stably stratified region of re- 
duced composition in light elements ? . Here we show that this layer can be generated by 
simultaneous crystallization and melting at the surface of the inner core, and that a transla- 
tional mode of thermal convection in the inner core can produce enough melting and crystal- 
lization on each hemisphere respectively for the dense layer to develop. The dynamical model 
we propose introduces a clear asymmetry between a melting and a crystallizing hemisphere 
which forms a basis for explaining the East- West asymmetry. The present translation rate is 
found to be typically 100 Ma for the inner core to be renewed entirely, which is one to two 
orders of magnitude faster than the growth rate of the inner core's radius. The resulting 
strong asymmetry of buoyancy flux caused by light elements is anticipated to have an impact 
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on the dynamics of the outer core and on the geodynamo. 

The original observation 7 of seismic P-wave velocities slower than the adiabatic PREM 
model in the lower outer core has since been confirmed and incorporated in ID global models 
AK135 7 and PREM2 ? . That discrepancy from the adiabatic profile could result from a wrong in- 
terpretation due to the nearby complex inner core, as sensitivity kernels have a width of several 
hundreds of km at body waves frequency 7 , or might also be attributed to floating crystals 7 ' 7 . Gub- 
bins et al. 7 show that this last explanation is not possible but that the observed seismic velocities 
can be explained by a stratification in light elements (and temperature). However, the stratifica- 
tion mechanism by crystallization and melting of crystals at different depths is not completely 
elucidated. 

We propose that a dense layer can develop when melting and crystallization only occur at the 
inner core boundary (ICB). Where crystallization takes place, light elements are released providing 
light fluid and where melting takes place, dense fluid is produced. It is possible to quantify these 
effects in terms of flux of buoyancy. Let us denote Ap that fraction of density jump across the 
ICB due to composition partition between solid and liquid phases. For a rate of crystallization 
V (respectively melting), the buoyancy flux is Ap g c V (respectively — Ap g c V), where g c is 
the magnitude of gravity 7 on the ICB. The idea is that part of the heavy fluid would remain at 
the bottom, while the rest would be entrained by the light fluid. Conversely, part of the light 
fluid would mix with the dense fluid in the dense layer while the rest would cross the dense layer 
and contribute to convection within the main part of the outer core. This idea has been validated 
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experimentally as follows. 

The experiments consist in injecting simultaneously a constant flux of light and dense fluid 
at the bottom of a fluid cavity. The cavity is a box of perspex 20 cm high and with a 15 cm x 15 
cm horizontal cross-section. It is initially filled with salted water (c wt % NaCl). At the bottom 
of the cavity, there is a porous layer (sponge) below which the cross-section is divided into two 
disconnected parts: on one side light fluid is injected (cj < c wt % NaCl) and on the other side 
heavy fluid is injected (c h > c wt % NaCl). Both density differences c — q and c h — c and 
both flow rates are controlled and set constant during the experiment. The injections of fluids start 
simultaneously through pipes from reservoirs with the desired concentration. The excess of fluid 
is removed through an overflow at the top of the cavity. 

The geophysically relevant case is when the positive buoyancy flux exceeds the negative 
one since the inner core is growing on average. When the negative buoyancy flux induced by the 
heavy fluid is less than 80 % in amplitude that of the light fluid, no dense layer is observed: the 
entrainment caused by the rise of light plumes is sufficient to mix the heavy fluid as it is released 
by the bottom boundary. However, when that heavy buoyancy flux is more than 80 % that of the 
light buoyancy flux, a growing dense layer forms at the bottom of the cavity. It has been observed 
experimentally that the condition for the existence of the dense layer is really a condition on the 
buoyancy fluxes, as described above, and not a condition on the volume flow rates nor on the 
density differences between the fluids. This constitutes a justification for the relevance of such a 
convection experiment as a model of a melting/crystallization process for the inner core. 
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On Fig. 1, an experimental run is shown. This experiment corresponds to a case where 
the heavy fluid buoyancy flux was 83 % that of the light fluid. The initial concentration and 
concentrations of the dense and light injected fluids were c = 4 wt %, c h = 6 wt % and q = 
1.65 wt % NaCl respectively. The volume flow rate of the dense fluid was Q h = 3.9 10~ 7 m 3 s _1 
and that of the light fluid Qi = 4.0 10 _7 m 3 s _1 . That experiment has been run twice under the same 
conditions: in the first instance, the injected dense fluid was colored with potassium permanganate 
and photographs of the setup have been taken at different times after the beginning of the injections. 
A dense colored layer forms at the bottom and its thickness grows linearly with time. It is also 
possible to see convection plumes going up on the right-hand side, carrying along some of the 
heavy colored fluid in the upper part of the cavity. In the second instance, the synthetic schlieren 
method has been used ? ? , providing a quantitative two-dimensional field of refraction index, hence 
visualizing concentration gradients: their horizontal components are shown on the middle row of 
Fig. 1 showing convection plumes of light fluid on the right-hand side of the cavity, while their 
vertical components are shown on the bottom row visualizing the dense layer and its growth. The 
concentration field is computed from its gradient, and averaged along the horizontal direction: the 
resulting stratification profile is shown on Fig. 2. There is clearly a region of stratified fluid, above 
which density is nearly uniform. The thickness of this layer is growing linearly with time, its 
volume being 50 % to 90 % that of the total volume generated by the light and heavy fluxes. 

Melting part of the inner core at a significant rate is difficult while it is crystallizing on 
average as a result of secular cooling. The most plausible way is that a topography is formed 
dynamically on the ICB so that the temperature of adjacent fluid of the outer core exceeds the 
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melting temperature. That excess temperature is then responsible for heat transfer from the outer 
core to the ICB, providing latent heat for fusion: hence topography can be related to the rate of 
melting. 



The dynamical model we put forward to account for significant melting on the ICB results 
from the combination of three physical elements: thermal state of a superadiabatic inner core, grav- 
itational equilibrium and finite heat exchange of latent heat with the outer core. In superadiabatic 
conditions, a uniform velocity in the inner core V, say from West to East along the rr-axis (see 
Fig. 3), generates a global superadiabatic temperature gradient in the same direction proportional 
to the residence time in the inner core, hence inversely proportional to V, and proportional to a 
positive source term S ~ 1CT 15 K s _1 defined from secular cooling and thermal conduction along 
the adiabat (see Methods Section and reference ? ): 

^ = ^ (1) 

where is the temperature relative to the adiabat T ad in the inner core anchored to the ICB ? . 
It follows from the volume expansion coefficient ? a = 1.1 10~ 5 K _1 and inner core density 
(on the ICB ? ) p s = 1.28 10 4 kgm~ 3 , that there exists a density gradient —ap s dQ/dx. The 
resulting gravity field and density distribution generate unbalanced forces on the inner core, so that 
it is displaced a distance 5 in the re-direction. In the Methods section, the gravitational field and 
potential associated with this mass distribution are derived, from which it is possible to calculate 
the net gravitational force F G exerted on the inner core and the net pressure force F P exerted by 
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the outer core on the inner core: 

f g + f p = —^-g p e c 3 



' dQ c 2 



(2) 



where Q is the universal gravitational constant, p e is the outer core density on the ICB ? and is 
the unit vector in the direction of the temperature gradient. The equilibrium condition that both 
forces balance provides the shift 5 as a function of the thermal gradient dQ/dx: 

5 = — ^Tl — (3) 
5 (Ps - Pi) 

Then, the displacement 5 is associated with a non uniform pressure distribution on the ICB (see 
Methods Section) hence to a small temperature departure 5T from the adiabat (see Fig. 4): 

5T = p e g c 5 cos 6 (m P - m ad ), (4) 

where m P = 8.5 1CT 9 K Pa -1 is the Clapeyron slope ? , m ad = (aT ad ) / (pc p ) = 6 1CT 9 K Pa -1 
is the adiabatic gradient and g c = Q^fp s c is gravity on the ICB. That departure is accomo- 
dated by a thermal boundary layer in the outer core, with a corresponding heat transfer of typi- 
cal magnitude u' c p 5T, where u' = 1CT 4 m s _1 is a typical velocity scale in the outer core and 
c p = 850 J kg -1 Kr 1 is the specific heat capacity ? . That heat transfer must be balanced by the 
release or absorption of latent heat: 

L V cos 9 = u' c p ST, (5) 

where L = 900 kJ kg -1 is the latent heat coefficient ? ? . Finally, combining equations ([!]), Q 
and (|5]), one can express the translational velocity: 

15 L (p s - p e ) 
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Depending on the heat flux at the CMB, the history of the inner core shows a first phase 
dominated by growth c ~ c _1 , followed by the development of the translational instability (see 
Supplementary information) when its radius was around 400 km, leading to the dominant present 
translation V ~ c 3 / 2 of order 5 1CT 10 m s _1 while the growth rate is of order 1CT 11 m s _1 (Fig. 5). 

The latter scaling law implies that the translational convection is faster along a long axis 
of the inner core oblate spheroid (see Supplementary information), i.e. perpendicular to the ro- 
tation axis. It follows that the temperature gradient is preferrentially aligned with such a long 
axis, which iagain reinforces convection in that direction. Moreover, the Earth's aspherical mass 
distribution which has essentially a degree 2, order 2 geometry ? is responsible for elongating 
slightly the inner core along an East- West axis and induces a degree 1 translational convection 
in the inner core through a bifurcation produced by instability (see supplementary information). 
We propose that the translational flow has a West to East orientation, hence being responsible for 
the observed hemispherical asymmetry of the inner core: grain growth during the transit from the 
western hemisphere to the eastern hemisphere may explain the difference in seismic properties ? . 
The temperature difference of a few K between both hemispheres is another source of asymmetry. 

According to our experiments, a melting rate above 80 % that of the crystallization rate is 
necessary for a dense layer to form, which geometrically implies that the translation velocity V is 
more than 20 times that of the inner core growth rate. From Fig. 5, this happens only when the 
CMB heat flux exceeds 10 TW, and only since the inner core radius was 1 100 km, some 200 Ma 
ago. Extrapolating from our experiments, evaluating 50 % of the volume of melt produced since 
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then would corresponds to a layer of thickness 250 km. The experimental excess concentration 
is found to be 10 % of the concentration difference between light and heavy injected fluids. In 
the Earth's core where concentration in light element is of order 10 %, a concentration difference 
around 1 % across the dense layer is expected. This is indeed coherent with the observed seismic 
velocities ? . 

Our convection mechanism ignores deformation in the inner core and compositional buoy- 
ancy. With a finite effective viscosity, temperature variations along gravity isopotentials induce an 
internal flow with deformation that affects the translational mode. We have estimated that an in- 
ternal flow is weak compared to translation for an effective viscosity above 10 18 Pa s. Enrichment 
in light elements of the outer core (a few percents) has been invoked ? ? to imagine a stabilizing 
mechanism for convection in the inner core. This is however very speculative, as the fraction of 
light elements incorporated in the inner core may have decreased more rapidly than the increase in 
outer core composition since gravity on the ICB is getting larger, reinforcing compaction. 

Invoking an excessively asymmetric buoyancy flux on the ICB deserves further study rela- 
tive to the dynamics of the outer core and the geodynamo. The stratified layer is expected to be 
dynamically isolated and to act as a filter between the inner core and the rest of the outer core, but 
there might subsist some hemispherical asymmetry in the outer core dynamics. 

METHODS SUMMARY 

The mode of convection associated to the translation of the inner core is not standard. Therefore, 
it is presented in the Methods section. Thermal buoyancy is the driving force, however, unlike 
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classical convection, the damping is not due to viscous and/or thermal diffusion. Damping is set 
by the capacity of the outer core to extract or supply latent heat on the ICB. 

Supplementary information is linked to the online version of the paper at www.nature.com/nature 
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Figure 1 Visualization of the growth of a dense layer in an experimental run, using (a) 
dye injection and (b) horizontal and (c) vertical density gradients. The experimental cavity 
is initially filled with a 4 wt% NaCI water solution. From t — 0, a constant flux of 1.65 
wt% NaCI solution is injected at the bottom on the right-hand side while a 6 wt% NaCI 
solution is injected on the left-hand side. The dense fluid is colored with permanganate 
potassium on the top row, visualizing a growing dense layer at the bottom, at four different 
times after the beginning of the injection. The synthetic schlieren method is used in a 
second identical experiment: the horizontal gradient of refraction index on the middle row 
highlights the convective plumes and the vertical gradient on the bottom row reveals the 
dense layer. 

Figure 2 Evolution of the concentration profile during the growth of a dense layer. The 
concentration field is extracted from the gradient of refraction index. It is averaged along 
the horizontal direction and shows the evolution of the dense layer. 

Figure 3 A schematic representation of the translational convective mode. The centre of 
the inner core O, is shifted by a distance 5 away from the centre of the Earth C, which 
would be its equilibrium position if its density were uniform. That shift causes thermody- 
namic disequilibrium at the ICB, generating melting on one side and crystallization on the 
other one. Hence a uniform flow exists in the inner core: in the case of a superadiabatic 
regime, a gradient of temperature develops as represented in grey scale. Its associated 
changes in density and gravitational potential lead to a new mechanical equilibrium for 
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the inner core, corresponding indeed to a shift in position in the same direction as initially 
assumed. 



Figure 4 Thermal departure from the adiabat due to the displacement of the inner core 
and heat transfer at the ICB. A thermal boundary layer forms in the outer core to adjust to 
the different radii of the ICB on the melting and crystallization sides. 



Figure 5 Growth rate of the radius of the inner core and uniform convective velocity as 
functions of the inner core radius. They are plotted for different values of the heat flux at 
the CMB. Light solid lines : mean solidifcation rate of the inner core. Light dash-dotted 
lines : translation velocity calculated with the assumption of a constant S. Heavy solid 
lines : translation velocity with S(t) calculated (see supplementary information) from the 
core thermal evolution model of ? . 
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